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1. Introduction 

The search for a quantum theory encompassing gravity has been a major issue in 
theoretical physics for the last 60 years. Nonetheless, until recently quantum gravity was 
largely relegated to the realm of speculation due to the complete lack of observational 
or experimental tests. In fact the traditional scale of quantum gravitational effects, the 
Planck scale Mpi = 1.2 x 10^^ GeV/c^, is completely out of reach for any experiment or 
observation currently at hand. This state of affairs had led the scientific community to 
adopt the "folklore" that testing quantum gravity was completely impossible. However 
the last decade has seen a dramatic change in this respect, and nowadays one can 
encounter a growing literature dealing with tests of possible predictions of various 
quantum gravity models [1, 2, 3, 4]. This field goes generically under the name of 
"quantum gravity phenomenology" . 

Among the several generic predictions associated with quantum gravity models, the 
possibility that Planck-scale physics might induce violations of Lorentz invariance has 
played a particularly important role [1, 2]. Generically any possible "discreteness" or 
"granularity" of spacetime at the Planck scale seems incompatible with strict Lorentz 
invariance (although some particular quantum gravity theories might still preserve it, 
see e.g. [5]) as larger and larger boosts expose shorter and shorter distances. Actually 
we now have a wealth of theoretical studies — for example in the context of string field 
theory [6, 7], spacetime foam scenarios [8], semiclassical calculations in loop quantum 
gravity [9, 10], DSR models [11, 12, 13], or non-commutative geometry [14, 15, 16, 17], 
just to cite a few — all leading to high energy violations of Lorentz invariance. 

Interestingly most investigations, even if they arise from quite different fundamental 
physics, seem to converge on the prediction that the breakdown of Lorentz invariance 
can generically become manifest in the form of modified dispersion relations exhibiting 
extra energy- dependent or momentum-dependent terms, apart from the usual quadratic 
one occurring in the Lorentz invariant dispersion relation: 

E^^m^ c'+p" c\ (1) 

In the absence of a definitive theory of quantum gravity it became common to adopt, 
in most of the literature seeking to put these predictions to observational test, a purely 
phenomenological approach, i.e., one that modifies the dispersion relation by adding 
some generic momentum- dependent (or energy-dependent) function F{p, c, Mpi) to be 
expanded in powers of the dimensionless quantity p/(Mpi c). Hence the ansatz reads: 

= ^2 ^4 ^ p2 ^2 ^ ^^^^ ^p^^ . ^2) 

oo 

= m2c"+p2c2 + J^t^^p'^; (3) 

n=l 

= + p" + c4^, M^,p/c + rj^p'/c' -h J] ^ 1 ; (4) 

I n>3 ^^''Pl J 

where all the zUn carry appropriate dimensions, and in contrast the ?7„ are chosen to be 
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dimensionless. Since these dispersion relations are not Lorentz invariant, it is necessary 
to specify the particular inertial frame in which they are given, and generally one chooses 
the CMB frame. Finally note that we have assumed rotational invariance and hence 
Lorentz violation only in the boost subgroup. This is motivated by the idea that Lorentz 
violation may arise in quantum gravity from the presence of a short distance cutoff. 
Moreover it is very difficult to conceive a framework where a breakdown of rotational 
invariance does not correspond to a violation of boost invariance as well. 

Of course merely specifying a set of dispersion relations is not always enough to 
place significant constraints — as most observations need at least some assumption 
on the dynamics for their interpretation. In fact most of the available constraints 
are extracted from some assumed model "test theory". Although several alternative 
scenarios have been considered in the literature, so far the most commonly explored 
avenue is an effective field theory (EFT) approach (see e.g., [18] for a review focussed on 
this framework, and [19, 20, 21, 22, 23, 24, 25] for some of the primary literature). The 
main reasons for this choice can be summarized in the fact that we are very familiar with 
this class of theories, and that it is a widely accepted idea (although not unanimously 
accepted, see e.g. [26]) that any quantum gravity scenario should admit a suitable EFT 
description at low energies, jj All in all, the standard model of particle physics and 
general relativity itself, (which arc presumably not fundamental theories) are EFTs, as 
are most models of condensed matter systems at appropriate length and energy scales. 
Even "fundamental" quantum gravity candidates such as string theory admit an EFT 
description at low energies (as perhaps most impressively verified in the calculations of 
black hole entropy and Hawking radiation rates). 

The EFT approach to the study of Lorentz violations has been remarkably useful 
in the last decade. Nowadays the best studied theories incorporating Lorentz violations 
are EFTs where Lorentz violations are associated either with renormalizable Lorentz- 
violating operators (mass dimension four or less), or sometimes with higher-order 
Lorentz-violating operators (mass dimension five and six or greater, corresponding to 
order and p'^ and higher deviations in the dispersion relation (4)). The first approach 
is generally known as the Standard Model Extension [27], while the second has been 
formalized by Myers and Pospclov [28] in the form of QED with dimension five Lorentz- 
violating operators (order p^ deviations in the dispersion relation of equation (4)). In 
both cases extremely accurate constraints have been obtained using a combination of 
experiments as well as observation (mainly in high energy astrophysics). See e.g., [1, 18]. 

In the present article we wish to focus on the non-renormalizable EFT with Lorentz 
violations developed in [28] , and subsequently studied by several authors. In spite of the 
remarkable success of this proposal as a "test theory" , it is interesting to note that there 
are still significant open issues concerning its theoretical foundations. In particular, let 
us now focus on two aspects of this approach that have spurred some debate among the 

II It is true that, e.g., non-commutative geometry can lead to EFTs with problematic IR/UV mixing, 
however this more likely indicates a physically unacceptable feature of such specific models, rather than 
a physical limitation of EFT. 
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quantum gravity phenomenology community. 

The naturalness problem: Looking at the dispersion relation (4) it might seem that the 
deviations linear and quadratic in p are not Planck suppressed, and hence are always 
dominant (and grossly incompatible with observations). However one might hope that 
there will be some other characteristic QFT mass scale -C Mpi (i.e., some particle 
physics mass scale) associated with the Lorentz symmetry breaking which might enter in 
the lowest order dimensionless coefficients 771^2, which will be then generically suppressed 
by appropriate ratios of this characteristic mass to the Planck mass. Following the 
observational leads one might then assume behaviours like 771 oc (/i/Mpi)°'+-^, and 
?72 oc (/i./Mpi)'^ where a > 1 is some positive power (often taken as one or two). 
Meanwhile no extra Planck suppression is assumed for the higher order rjn coefficients 
which are then naturally of order one. Note that such an ansatz assures that the Lorentz 
violation term linear in the particle momentum in equation (4) is always subdominant 
with respect to the quadratic one, and that the Lorentz violating term cubic in the 
momentum is the less suppressed of the higher order ones. ^ If this is the case one 
will have two distinct regimes: For low momenta p/{M-p\c) <^ {ji/MpiY the lower-order 
(quadratic in the momentum) deviations in (4) will dominate over the higher-order 
(cubic and higher) ones, while at high energies p/{M-pic) ^ (///Mpi)'^ the higher order 
terms (cubic and above in the momentum) will be dominant. 

The naturalness problem arises because such a line of reasoning does not seem to 
be well justified within an EFT framework. In fact we implicitly assumed that there 
are no extra Planck suppressions hidden in the dimensionless coefficients 7]^ with n > 3. 
Indeed we cannot justify why only the dimensionless coefficients of the n < 2 terms 
should be suppressed by powers of the small ratio /i/Mpi. Even worse, renormalization 
group arguments seem to imply that a similar mass ratio, /i/Mpi would implicitly be 
present in all the dimensionless n > 3 coefficients — hence suppressing them even 
further, to the point of complete undetectability. Furthermore it is easy to show [29] 
that, without some protecting symmetry, it is generic that radiative corrections due 
to particle interactions in an EFT with only Lorentz violations of order n > 3 in (4) 
for the free particles, will generate n = 1 and n = 2 Lorentz violating terms in the 
dispersion relation which will then be dominant. 

The universality issue: The second point is not so much a problem, as an issue of debate 
as to the best strategy to adopt. In dealing with situations with multiple particles one 

^ Of course this is only true provided there is no symmetry like parity that automatically cancels all 
the terms in odd powers of the momentum. In that case the least suppressed Lorentz-violating term 
would be that quartic in the momentum. 

+ A symmetry which could play a protective role for the lowest-order operators as indeed been 
suggested. In [30] it was shown that the dual requirements of supersymmetry and gauge invariance 
permit one to add to the SUSY standard model only those operators corresponding to n > 3 terms 
in the dispersion relation. However it should be noted that in [30] the 773 coefficient carries a further 
suppression of order nt^/M^ when compared to (4). 
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has to choose between the case of universal (particle- independent) Lorentz violating 
coefficients ry^, or instead go for a more general ansatz and allow for particle-dependent 
coefficients; hence allowing different magnitudes of Lorentz symmetry violation for 
different particles even when considering the same order terms (same n) in the 
momentum expansion. The two choices are equally represented in the extant literature 
(see e.g. [31] and [2] for the two alternative ansatze), but it would be interesting to 
understand how generic this universality might be, and what sort of processes might 
induce non-universal Lorentz violation for different particles. 



To shed some hght on these issues it would definitely be useful to have something 
that can play the role of test-bed for some of the ideas related to the emergence and 
form of the modified dispersion relations of equation (4). In this regard, herein we 
will consider an analogue model of emergent spacetime, that is, a condensed matter 
system which admits excitations whose propagation mimics that of quantum fields on a 
curved spacetime [32].* Indeed it is well known that the discreteness at small scales of 
such systems shows up exactly via modified dispersion relations of the kind described 
by equation (4), and one may hope that the complete control over the microscopic 
(trans-Planckian) physics in these systems would help understanding the nature of the 
issues discussed above [36, 37]. For example, we remind the reader that the Bogolubov 
quasi-particle spectrum for excitations of a Bose-Einstein condensate 

uj''^clk'' + clk^/K\ (5) 

where Cg is the speed of sound and K is determined by the effective Compton 
wavelength [36, 37] 

27r „ mCs 

K = = 27r— ^, (6) 

'^Compton I'' 

i.e., it is set by the the mass of the fundamental bosons forming the BEC. Note also 
that we do not encounter odd powers of the momentum in (5), as the system is by 
construction invariant under parity. 

This system (which has been extensively studied as an analogue model of 
gravity [36, 37, 38], in particular with reference to the simulation of black holes via 
supersonic flows), provides a simple and explicit example of the high-energy breakdown 
of "Lorentz invariance" [36, 37] with a dispersion relation of the form (4) which 
interpolates between a low-energy "massless" relativistic regime 

^ cie, {k<^K)- (7) 

and a non-relativistic, approximately Newtonian, high-energy regime 

ujp^Csk'^/K, (k:>K). (8) 

* For other rather distinct views on "emergent spacetime" , with rather different aims, see for instance 
the articles by C D Proggatt and H B Nielsen [33], by J D Bjorken [34], and by R B Laughlin [35]. 
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Here K is hence the scale of the violation of the Lorentz invariance and as such in the 
language of quantum gravity phenomenology it is the analogue of the Planck scale, ji 

Unfortunately it is easy to realize that this particular system is too simple in order 
to mimic the salient issues of quantum gravity phenomenology. Actually all of the 
analogue models currently at hand have a problem in this sense, either because they do 
not provide the same dispersion relations for different excitations even at low frequencies 
(in the words of [39, 40], "no mono-metricity" ) , or because they deal with just one single 
kind of excitation (hke the 1-BEC model just discussed) — in which case it is impossible 
to say anything about either naturalness or universality. [Regarding the naturalness 
problem, if there is only one type of excitation present in the system then one cannot 
identify a 772 modification at the quadratic (p^) order, given that it shows up only via 
differences in the limit speed when comparing different particles.] 

However it was recently realized that experimentally available systems of two 
coupled BECs are much richer in their spectrum, and allow the simulation of an analogue 
spacetime where two different particles coexist and interact through mode-mixing. In 
particular in [39, 40] it was shown that for a 2-BEC system there are ways of modifying 
the excitation spectrum (5) in order to add a "mass" term, and the analysis of the 
present article will build on those two articles. (For a somewhat related though distinct 
application of 2-BEC systems to analogue models, see [41].) 

We shall consider the special case of a homogeneous two-component BEC subject to 
laser-induced coupling. This system exhibits a rich spectrum of excitations, which can 
be viewed as two interacting phonon modes (two quasiparticlc modes). We study the 
conditions required for these two phonon modes to share the same "special relativity" 
metric in the hydrodynamic limit (effectively the low-energy limit), and find that (in this 
limit) the two phonons respectively exhibit a massive and massless Lorentz invariant 
dispersion relation. We then consider the high energy limit of the system, that is, 
situations where the "quantum pressure" term, which in a 1-BEC system is is at the 
origin of the k'^ Lorentz violation in (5), can no longer be neglected. 

Though much of the underlying physics is similar to that of [39, 40], the central 
thrust of the argument is different. In those articles one was always working in the 
hydrodynamic limit, often with inhomogeneous backgrounds, and seeking to extract 
a curved spacetime metric. In the current article we are working with homogeneous 
backgrounds, staying as close as possible to "special relativity" , and specifically probing 
the possible breakdown of Lorentz invariance by going beyond the hydrodynamic 
approximation. Thus many of the issues that normally are central to the discussion 
of analogue models (such as the existence of a curved effective spacetime, analogue 
horizons, causal structure, analogue Hawking radiation, the simulation of cosmological 

Let us stress that while it is a standard assumption in quantum gravity phenomenology to identify 
the scale of Lorentz violation with the Planck scale, it is not a priori necessary that the two must 
exactly coincide. In the discussion below, what we shall call the effective or analogue "Planck scale" 
has always to be interpreted as the scale of Lorentz breaking; given that this is the only relevant scale 
in our discussion. 
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spacetimes, and the like [42, 43, 44]), here are of at best peripheral interest. The 
analogue spacetimes we are interested in are all flat. Finally we also stress that our main 
thrust is here to (eventually) learn lessons about how real physical quantum gravity 
might work — we are not particularly concerned about the experimental laboratory 
realizability of our specific analogue system. (Readers more interested in specific 
condensed-matter aspects of 2-BEC systems and their excitations might consult, for 
instance, references [45, 46, 47], and references therein.) 

The rest of the paper will be organized as follows. In section 2 we shall describe 
the general equations for the propagation of excitations of the 2-BEC background, while 
section 3 then discusses the notion of "healing length" in a 2-component BEC. Then in 
section 4 we shall explicitly consider the low energy limit when the quantum potential 
{i.e., the ultraviolet physics due to the atomic nature of the condensate) is negligible. 
This limit will allow us to identify for which combination of the microscopic parameters 
this system is indeed an analogue model of gravity, i.e., it is characterized by a single 
metric for all the excitations. In section 5 we then move on to explore how UV effects, 
embodied by the quantum potential, introduce Lorentz violations, and in section 6 we 
study the implications for quantum gravity phenomenology. In the final discussion, 
section 7, we consider the lessons one can draw from this analogue model about how 
low energies Lorentz violations can be protected in realistic situations. 



2. Phonons in two-component BECs 



Two BECs interacting with each other, and coupled by a laser-driven coupling, can 
usefully be described by the pair of Gross-Pitaevskii equations [39, 40]: 



ihdt^A 



2 niA 
2mB 



+ - /X + Uaa 
V^ + Vb-h + Ubb 





^ + Uab 






1' + Uab 





(9) 
(10) 



We permit ruA ^ tub va. the interests of generality, (although ruA ~ raB in all 
currently realizable experimental systems) , and note that A can take either sign without 
restriction. In the specific idealized case A = and Uab = we have two inter- 
penetrating but non-coupled BECs, each of which separately exhibits a Bogolubov 
spectrum with distinct values of the Lorentz breaking scale K (since uia 7^ w^b). Once 
the BECs interact the spectrum becomes more complicated, but that is exactly the case 
we are interested in for this article. 

To analyze the excitation spectrum we linearize around some background using: 

= ^Pxo + epxi e'(^^°+"^^i) for X = A,B. (11) 

Because we are primarily interested in looking at deviations from special relativity (SR), 
we take our background to be homogeneous [position-independent], time-independent, 
and at rest [vaq = = vbo]- That is, we will be dealing with an analogue of 
flat Minkowski spaccetime. We also set the background phases equal to each other. 
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^Ao = ^Bo- This greatly simplifies the technical computations. ff Then the linearized 
Gross-Pitaevskii equations imply: 

Oai = -—^ PAi-^^^^ PB1 + -^^Qai{pai), (12) 

-Qm(pm), (13) 



Uaa 


Uab 


Pbi + 


h 


h 


h 


1 


Ubb 


Uab 


Pai + 


h 




h 


2mB 



and 



TBI 



h 2A 

Pai = Pao'^^Oai + -^VPao Pbo {Obi -Oai), (14) 

iJiA ri 

h 2A 

Pbi = Pbo^^Obi + -^^fpmPm {Oai - Obi)- (15) 

TUB ri 



Here we have defined 



~ A TT ^VPaoPbo 1 
Uaa — Uaa — — . 3 = Uaa —■, 



2 ^Pao 2 p^Q 

A v^Pao _ jj ^\/pAQpm 1 

2v^ " ^ ;^o' 



fr —TT ^ 7-r ^^y/HAOHBO -L . . 

t^BB — t^BB — — , 3 — t^BB 2~' v 



~ _ A 1 _ AVPaoPbo 1 

(7ab — t/AB + — : — Uab + 



2 ^/pAo'pBo 2 p AO Pbo 

and furthermore we define Qxi as the second-order differential operator obtained from 
linearizing the quantum potential: 

'^^^^^ 2mx V y 2mx V VPxo + epx, J ^ ' 



2 



(Qxo(pxo) + e Qxi{pxi)) ■ (18) 

The quantity Qxo{pxo) corresponds to the background value of the quantum pressure, 
and contributes only to the background equations of motion — it does not affect the 
fiuctuations. Now in a general background 

9x.fe.) - ^ I - <^'>'^°)>'^° - YJ^V + p,,. (19) 

2 I Pxo Pxo Pxo J 

Given the homogeneity of the background appropriate for the current discussion this 
simplifies to 

Qxi{pxi) = TT^VVxi. (20) 
2pxo 

The set of first-order partial differential equations relating the phase fiuctuations 
and density fiuctuations can be written in a more concise matrix form. First let us 
define a set of 2 x 2 matrices, starting with the coupling matrix 

E^E + X, (21) 
ft For some of the additional complications when background phases are unequal see [39]. 



Analogue quantum gravity phenomenology from a two-component BEC 
where 



„ _ 1 



Uaa Uab 
Uab Ubb 



and 



Qai 
niA 







Qbi 



h 
i 



rUA PAO 







1 



rriB pBO 



V' = -X v 



A second coupling matrix can be introduced as 
2A y/PAo Pbo 



h 



1 -1 
■1 +1 



It is also useful to introduce the mass-density matrix D 

^ _ ^ pAo/rriA 

PBo/ruB 

Now define two column vectors 9 — [9ai, Obi]^ and p — [pai, Pbi]^- 



(22) 



(23) 



(24) 



(25) 



Collecting terms into a 2 x 2 matrix equation, the equations for the phases (12) 
and densities (14) become 

0^-Ep, (26) 



'p= -D v^e + . 



(27) 



Equation (26) can now be used to eliminate p in equation (27), leaving us with a single 
matrix equation for the perturbed phases: 

dt (e-^ ^DW^e-Ae (28) 

This is an integro-differential equation [since S is a matrix of differential operators] which 
is a second-order differential equation in time, but an integral equation (equivalently, 
an infinite-order differential equation) in space. 

We now formally construct the operators S^/^ and S~^/^ and use them to define a 
new set of variables 

e = S-^/' e, (29) 

in terms of which the wave equation becomes 

d^e = [^v^ - A] E'/^'je, (30) 

or more explicitly 

d^e = {[s - xvT^' [Dv^ - A] [2 - xvT^'} &■ (si) 

This is now a (relatively) simple PDE to analyze. Note that whereas the objects S^/^ 
and S~^/^ are 2x2 matrices whose elements are pseudo-differential operators, as a 
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practical matter we never have to descend to this level of technicality. Indeed, it is 
computationally efficient to directly go to the eikonal limit where 

S^S + XA;^ (32) 

This finally leads to a dispersion relation of the form 

det{o^' I - [S + Xk^Y'^ [De + A] [S + X^Y'^] = , (33) 

and "all" wc need to do for the purposes of this article, is to understand this quasiparticle 
excitation spectrum in detail. 

3. Healing length 

Note that the differential operator Qxi that underlies the origin of the X contribution 
above is obtained by linearizing the quantum potential 

which appears in the Hamilton-Jacobi equation of the BEC flow. This quantum 
potential term is suppressed by the smallness of h, the comparative largeness of nix, and 
for sufficiently uniform density profiles. But of course in any real system the density of a 
BEC must go to zero at the boundaries of its EM trap (given that px = I'ipxix, t)P). In a 
1-component BEC the healing length characterizes the minimal distance over which the 
order parameter goes from zero to its bulk value. If the condensate density grows from 
zero to po within a distance ^ the quantum potential term (non local) and the interaction 
energy (local) are respectively £^kinetic ~ ^^/(2m^^) and ^'interaction ~ ATrh^apo/m. These 
two terms are comparable when 

e = (STTpoa)-'/', (35) 
where a is the s-wave scattering length defined as 

Note that what we call Uq in the above expression is just the coefficient of the non- 
linear self-coupling term in the Gross-Pitaevskii equation, i.e., just Uaa or Ubb ff we 
completely decouple the 2 BECs {Uab = A = 0). 

From the definition of the healing length it is hence clear that only for excitations 
with wavelengths much larger than the healing length is the effect of the quantum 
potential negligible. This is called the hydrodynamic hmit because the single-BEC 
dynamics is described by the continuity and Hamilton-Jacobi equations of a super- 
fiuid, and its excitations behave like massless phononic modes. In the case of excitations 
with wavelengths comparable with the healing length this approximation is no longer 
appropriate and deviations from phononic behaviour will arise. 

Such a simple discrimination between different regimes is lost once one considers 
a system formed by two coupled Bose-Einstein condensates. In fact in this case one is 
forced to introduce a generalization of the healing ^ length in the form of a "healing 
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matrix". Let us elaborate on this point. If we apply the same reasoning used above 
for the definition of the "healing length" to the 2-component BEC system (9, 10) we 
again find a functional form like that of equation (35) however we now have the crucial 
difference that both the density and the scattering length are replaced by matrices. In 
particular we generalize the scattering length a to the matrix A: 



A 



/niA 








Uaa 
Uab 



Uab 
Ubb 



'niA 






^/rriB 



Furthermore, from (35) a healing length matrix Y can be defined by 





Y 



-2 



2 



y/PAOrflA 

y/PBoiriB 



Uaa 
Uab 



Uab 
Ubb 



^ypAorriA 
y/PBorriB 



That is, in terms of the matrices we have so far defined: 



(37) 



(38) 



(39) 



We can now define "effective" scattering lengths and heahng lengths for the 2-BEC 
condensate as 



1 



and 



That is 



seff 



tr[^] 



rriAUAA + rnsUBB 



tr[y2] = tr[XS-i] 



h^[UBB/{mAPAo) + UAA/{mBPBo)] 



^{UaaUbb - U 



AB) 



^^[mAPAoUAA + msPBoUBB] 



(40) 



(41) 



(42) 



"imAmBpAopBo {UaaUbb - U\j^) 
Note that if the two components are decoupled and tuned to be equivalent to each other, 
then these effective scattering and healing lengths reduce to the standard one-component 
results. We shall soon see that sometimes it is convenient to deal with explicit formulae 
involving the "low level" fundamental quantities such as m^, m^, Uxyi etc., while often 
it is more convenient to deal with "high level" quantities such as ^eflf- 



4. Hydrodynamic approximation 

We are now interested in investigating the most general conditions (in the hydrodynamic 
limit) under which this two-EEC system can describe two phononic modes propagating 
over the same metric structure. The hydrodynamic hmit is equivalent to formally setting 
X — > so that S — > S. (That is, one is formally setting the healing length matrix to 
zero: y ^ 0. More precisely, all components of the healing length matrix are assumed 
small compared to other length scales in the problem.) 
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^.i. Fresnel equation 

The PDE (28) now takes the simphfied form: 

d^e = +[ED]- V¥ - [S A] ■ 9. (43) 

Since this is second-order in both time and space derivatives, we now have at least 
the possibihty of obtaining an exact "Lorentz invariance". We can now define the 
transformed variables 

and the matrices 

Q2 = A sV2; ^ hV2 d sV2; (45) 

so that 

d^e = +c^v^9 - n^e. (46) 

Then in momentum space 

cu^e = {C^ + 9 = H{k^) 9, (47) 
leading to the Presnel equation 

det{a;^ I - //(A;^)} = 0. (48) 



That is 



whence 



cu^ - u;^ tr[H{k^)] + det[H{e)] = 0, (49) 



2 _ tr[//(A;2)] ± ytr[//(F)]2 - 4 det[H{k^)] 



= ^ ^ " ^ ' (50) 

Note that the matrices Q^, Cq, and H{k^) have now carefully been arranged to be 
symmetric. This greatly simplifies the subsequent matrix algebra. Also note that the 
matrix H{k'^) is a function of /c^; this will forbid the appearance of odd powers of k 
in the dispersion relation — as should be expected due to the parity invariance of the 
system. 

4.2. Masses 

We read off the "masses" by looking at the special case of space-independent oscillations 
for which 

d^9 = -n^ 9, (51) 

allowing us to identify the "mass" (more precisely, the natural oscillation frequency) as 
"masses" oc eigenvalues of (E}^^ A E^^^) — eigenvalues of (S A). (52) 
Since A is a singular 2x2 matrix this automatically implies 

(jjj = 0; ujjj = ti{E A). (53) 
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So we see that one mode will be a massless phonon while the other will have a non zero 
mass. Explicitly, in terms of the elements of the underlying matrices 

u^h--^-^^^W^{UAA + UBB-2U^s} (54) 
nr 

so that (before any fine-tuning or decoupling) 
^11 = — { Uaa + Ubb - 2Uab 



Pao , Pbo 



It is easy to check that this quantity really does have the physical dimensions of a 
frequency. 

J!i..3. Conditions for mono-metricity 

We now want our system to be a perfect analogue of special relativity. That is: 

• We want each mode to have a quadratic dispersion relation; 

• We want each dispersion relation to have the same asymptotic slope. 

In order to find under which conditions these requirements are satisfied we will adopt 
the following strategy: Let us start by noticing that the dispersion relation (50) is of 
the form 

cu'^ — [quadratic^] ± ^/[quartic]. (56) 
The first condition imphes that the quartic must be a perfect square 

[quartic] = [quadratic2]^, (57) 

but then the second condition implies that the slope of this quadratic must be zero. 
That is 

[quadratic2](/c^) = [quadratic2](0), (58) 

and so 

[quartic] (A;^) = [quartic] (0) (59) 

must be constant independent of k"^, so that the dispersion relation is of the form 

o;^ = [quadratiCj](A;^) ± [quadratic2](0). (60) 

Note that this has the required form (two hyperbolae with the same asymptotes, and 
possibly different intercepts). Now let us implement this directly in terms of the matrices 
and M2. 

Step 1: Using the results of the appendix, specifically equation (A. 6): 

dei[H'^{k)] = det[Q^ + k^] (61) 

= det[Q'] - tr {n^ C^} e + det[Co'] {k^f. (62) 

(This holds for any linear combination of 2 x 2 matrices. Note that we apply trace 
reversal to the squared matrix Cq, we do not trace reverse and then square.) Since in 
particular det[Q^] = 0, we have: 

dei[H^{k)] = -tr d] k^ + det[C2] {k^f. (63) 
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Step 2: Now consider the discriminant (the quartic) 

quartic = ir[H{k'^)f - 4 det[H{k'^)] (64) 

= (tr[Q2] + ir[Cl] - 4 [-tr {Q^ C^} ^ + det[C^] (A;')'] (65) 
= tr[l]Y + {2tr[l]2]tr[Co'] + 4tr {fi^ 0^}}^ 

+ {tr[CoT-4det[Co^]}(fc^)2 (66) 
= tT[Q^]^ + 2{2tr {17" Co'} - tT[Q^]tT[C^]}k^ 

+ {tr[CoT-4det[Co']}(A;2)l (67) 

So in the end the two conditions above for mono-metricity take the form 

mono-metncity ^ <( _ ^^^^.j ^^^^.^ ^ ^ (68) 



Once these two conditions are satisfied the dispersion relation is 

2 tr[H(e)] ± tr[fi2] tr[Q2] ± tr[Q2] + tr[Co'] ^ 
^ = 2 = 2 

whence 

as required. One mode is massless, one massive with exactly the "mass" previously 
deduced. One can now define the quantity 

mn = lkoij/cl, (71) 

which really does have the physical dimensions of a mass. 

4-4- Interpretation of the mono-metricity conditions 

But now we have to analyse the two simplification conditions 

CI : tr[CoT - 4 det[Co'] = 0; (72) 
C2 : 2 tr c^j _ tr[Q^]tT[C^] = 0; (73) 

to see what they tell us. The first of these conditions is equivalent to the statement 
that the 2x2 matrix Cq has two identical eigenvalues. But since Cq is symmetric this 
then implies Cq = Cq I, in which case the second condition is automatically satisfied. (In 
contrast, condition C2 does not automatically imply condition CI.) Indeed if Cq = Cq I, 
then it is easy to see that (in order to make Cq diagonal) Uab — 0, and furthermore 
that 

Uaa Pao 2 Ubb Pbo 

= Co = . 74) 

niA TUB 

Note that we can now solve for A to get 

A = -l^pAo Pbo Uab, (75) 
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whence 



2 Uaa Pao + Uab Pbo Ubb Pbo + Uab Pao frrc\ 

Co = = , (76) 

ruA iriB 



and 

An An n r,nT I A 1-, I r I n Ar, I rin^ 1^1 

(77) 



2 ^PAOPBoUaB , jj . jj r,TT , TT 

uJn = ^ \ Uaa + Ubb - ^Uab + Uab 



Pao _^ / Pbo 



Pbo V Pao 



Note that (77) is equivalent to (55) with (75) enforced. But this then imphes 

^PAoPBoUaB ] jj , rr . tt P^ _^ PBO 

Pbo Pao 



^11 = ^2 1 + + ^^-B 



(78) 



Interpretation: Condition C2 forces the two low-momentum "propagation speeds" to 
be the same, that is, it forces the two O(fc^) coefficients to be equal. Condition CI is the 
stronger statement that there is no 0(/c^) (or higher order) distortion to the relativistic 
dispersion relation. 

5. Beyond the hydrodynamical approximation 

At this point we want to consider the deviations from the previous analogue for special 
relativity. To do so we have to reintroduce the linearized quantum potential Qxi in our 
equation (28) given that, at high frequencies, it is this term that is inducing deviations 
from the superfluid regime for the BEC system. 

5.1. Fresnel equation 

Our starting point is equation (33) which we rewrite here for convenience: 

u^e = |V5 + A: A;2 [D e + A] VE + X k^^ 9 = H{k^) 9. (79) 

This leads to the Presnel equation 

det{u;2 I - H{e)} = 0. (80) 

That is 

ou^ - tY[H{k^)] + det[H{e)] = 0, (81) 

whence 



ir[H{k^)] ± ^tr[H{k^)f - 4 det[H{k^\ 
^ = 2 ' ^ ' 

which is now of the form 

— [quartiC]^] ± [octic] . (83) 
And we can now proceed with the same sort of analysis as in the hydrodynamical case. 
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5.2. Masses 

The "masses" , defined as the zero momentum oscillation frequencies, are again easy to 
identify. Just note that 

H{k^ ^ 0) = n^, (84) 

and using the fact that dct(Q^) = one again obtains 

to^k^O) = {0,tr[n^]}. (85) 

So even taking the quantum potential into account we completely recover our previous 
results (53), (54), and (55). Of course this could have been predicted in advance by 
just noticing that the /c-independent term in the Fresnel equation is exactly the same 
mass matrix fi^ = S^/^ A S^/^ as was present in the hydrodynamical hmit. (That is, 
the quantum potential term X does not influence the masses.) 



5.3. Dispersion relations 

Let us start again from the general result we obtained for the dispersion relation (82). 
Differently from the previous case, when the hydrodynamic approximation held, we now 
have that the discriminant of (82) generically can be an eighth-order polynomial in k. In 
this case we cannot hope to recover an exact analogue of special relativity, but instead 
can at best hope to obtain dispersion relations with vanishing or suppressed deviations 
from special relativity at low k; possibly with large deviations from special relativity 
at high momenta. From the form of our equation it is clear that the Lorentz violation 
suppression should be somehow associated with the masses of the atoms niA/B- Indeed 
we will use the underlying atomic masses to define our "Lorentz breaking scale" , which 
we shall then assume can be identified with the "quantum gravity scale". The exact 
form and relative strengths of the higher-order terms will be controlled by tuning the 
2-BEC system and will eventually decide the manifestation (or not) of the naturalness 
problem and of the universality issue. 

Our approach will again consist of considering derivatives of (82) in growing even 
powers of k"^ (recall that odd powers of k are excluded by the parity invariance of the 
system) and then setting A; ^ 0. We shall compute only the coefficients up to order 
k^ as by simple dimensional arguments one can expect any higher order term will be 
further suppressed with respect to the k'^ one. 

We can greatly simplify our calculations if before performing our analysis we 
rearrange our problem in the following way. First of all note that by the cyclic properties 
of trace 

tr[//(A;2)] = tr[(L>A;2 + A) (S + k'^X)] (86) 
= tr[AS + P(i:>S + AX) + (A;2)2M] (87) 
= tr[sV2ASV2 + A;2(S1/2L>sV2 + X^^A^V^) + {efx^'''DX^/\ (88) 

We can now define symmetric matrices 

= E'/^AE'/^; (89) 
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With all these definitions we can then write 



ir[H{k'^)\ = tr 



where everything has been done inside the trace. If we now define 

then Hs{k^) is by definition both polynomial and symmetric and satisfies 

tr[H{e)]=tT[H,{k% 
while in contrast, 

det[ii'(A;2)] ^ det[Hs{k% 



But then 



Whence 



and at A; = 



dp 



da;2 



tr[Hs{e)] ± Vtr[//,(A;2)]2 - 4det[//(A;2)]J . 

tr[Hs{k^)]tT[H',{e)] - 2det'[H{e)] 



Vtr[i/,(A;2)]2 - 4det[ii"(A;2 



But now let us consider 



tr[C^] ± 



^,^tr[n^]tr[C']-2det'[H{k')]k^o 
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(90) 

(91) 

(92) 
(93) 

(94) 

(95) 
(96) 
(97) 

(98) 
(99) 



det[H{k^)] = det[{De + A) (S + k'^X)] 
^det[Dk^ + A] det[E + k'^X] 
= det[E'/\Dk^ + A)SV2] det[/ + k^E-'/^XE-'/^] 

where we have repeatedly used properties of the determinant. Furthermore 

det[/ + k^E-'/'XE-'/''] = det[/ + eE-'X] 

= det[I + ex'/'EX'/'] 
= det[/ + A;2f72], 

so that we have 

det[ii'(A;')] = det[n' + C^k^] det[I + k'^Y^ /2]. 

Note the the matrix is the "heahng length matrix" we had previously defined, and 
that the net result of this analysis is that the full determinant is the product of the 



(100) 
(101) 
(102) 

(103) 
(104) 
(105) 

(106) 
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determinant previously found in the hydrodynamic limit with a factor that depends on 
the product of wavenumber and healing length. 

But now, given our formula (A. 6) for the determinant, we see 

det'[//(A;2)] = (-tr(Q2c'2) + 2k^det[Cl\) det[7 + k''Y^l2] 

+ det[Q' + Cie] (-tr[y2] + A;Met[y'])/2, 



whence 



and so 



That is: 



da;2 



dA;2 
do;2 



dfc2 



= -iiin'^Cl), 

tr[f]2]tr[C2] + 2tr(^]2c'2^ 



tr[C^] ± ^tn^^i, 



tr[C2] ± <^ tr[C2] + 2 



tr[02 



21 , M^'ci) 



tr[r]2 



(107) 

(108) 

(109) 
(110) 



Note that all the relevant matrices have been carefully symmetrized. Also note the 
important distinction between Cq and C^. Now define 



then 



with 



da;2 



dA;2 



V2 



tr[C^]tr[r]^] + 2tr{Q^Ci) 
tr[fi2]tr[C2] 

Similarly, consider the second derivative: 

1 



tr(O^Cg 



(111) 



(112) 



(113) 



d(F)2 2 
± 



tr[i/;(A;2)] 

tr[g,(A;^)]tr[g;(P)] + tr[ff^(A:^)]tr[i/^(P)] - 2det^^[g(A;^)] 

v/tr[iJ,(A;2)]2 - 4det[if(A;2)] 

(tr[g,(A;^)]tr[i7^(A;^)] - 2dct^[iJ(P)])^ ' 
(tr[i/,(A;2)]2 -4det[i/(A;2)])3/2 



(114) 



whence 



d(fc2)2 



tr[Z^] ± 



tr[n^]tT[Z^] + tr[C^]^ - 2det^^[ij'(fc^)]fc^o 
tr[l]2] 



(tr[(]2]tr[C2] - 2det'[i7(/t2)]fc^o)'' 



tr[Q 



213 



(115) 
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The last term above can be related to da;^/dA;^, while the determinant piece is evaluated 
using 

2^^2 /ol (llg) 



d.ei"[H{e)] = (2det[Co']) det[/ + 



Therefore 



dei"[H{k'')] 



That is, (recalling tr[A] 



+ (-tr(Q2c'2) + 2A;Met[Co']) (-tr[y'] + k''det[Y^])/2 
+ det[Q' + Co'A;'] {det[Y^]/2) 

+ (-tr(Q'C'o') + 2A;Met[Co']) (-trfy^] + A;Met[y'])/2. 

(2det[Co^]) 

+ (-11(^2^2)) (-tr[F2])/2 + det[Q^] (det[F2])/2 

+ {-tT{n'C',)) (-tr[F^])/2. (117) 

tT[A]), 



det"[H{k')]k^o = (2det[C2]) - {tr{n'C',)) {tr[Y']), 



or 



Now assembling all the pieces 



(118) 



(119) 



d(F)2 



That is 

d^UJ^ 



tT[Z] ± 



tr[n2]tr[Z2] + tr[C2]2 + 2tr[CoC'2] ^ 2{tr{n^C^)) {tr[Y^]) 



tr[n2 



trp 



■213 



(120) 



2^2 



d{k 



and so 



fc^O 



tr[Z] ± 
^ tr[Q2] 



tr[l]2]tr[Z2] + tr[C2]2 + 2tr[C2q] + 2{tT{Q^C^)) (trfF^]) 



tr[n2 



(121) 



d(A; 



2^2 



1 
2 



tr[Z2] ± tr[Z2] ± 2 



tr[fi^, 



°Jtr[F2]± 



tr[C2]2-4det[C, 



tr[02 



tr[C- 



■212 



tr[l]2 



(122) 



With the above formula we have completed our derivation of the lowest-order terms 
of the generic dispersion relation of a coupled 2-BEC system — including the terms 
introduced by the quantum potential at high wavenumber — up to terms of order A;^. 
Prom the above formula it is clear that we do not generically have Lorentz invariance 
in this system: Lorentz violations arise both due to mode-mixing interactions (an effect 
which can persist in the hydrodynamic limit where Z ^ and Y ^ 0) and to the 
presence of the quantum potential (signalled by Z 7^ and Y ^ 0). While the mode- 
mixing effects are relevant at all energies the latter effect characterizes the discrete 
structure of the effective spacetime at high energies. It is in this sense that the quantum 
potential determines the analogue of quantum gravity effects in our 2-BEC system. 



Analogue quantum gravity phenomenology from a two-component BEC 



20 



6. The relevance for quantum gravity phenomenology 

Following this physical insight we can now easily identify a regime that is potentially 
relevant for simulating the typical ansatze of quantum gravity phenomenology. We 
demand that any violation of Lorentz invariance present should be due to the microscopic 
structure of the effective spacetime. This implies that one has to tune the system in 
order to cancel exactly all those violations of Lorentz invariance which are solely due to 
mode- mixing interactions in the hydrodynamic limit. 

Wc basically follow the guiding idea that a good analogue of quantum-gravity- 
induced Lorentz violations should be characterized only by the ultraviolet physics of the 
effective spacetime. In the system at hand the ultraviolet physics is indeed characterized 
by the quantum potential, whereas possible violations of the Lorentz invariance in the 
hydro dynamical limit are low energy effects, even though they have their origin in the 
microscopic interactions. We therefore start by investigating the scenario in which the 
system is tuned in such a way that no violations of Lorentz invariance are present in the 
hydrodynamic hmit. This leads us to again enforce the conditions CI and C2 which 
corresponded to "mono-metricity" in the hydrodynamic limit. 

In this case (110) and (122) take respectively the form 



da;2 



and 



dA; 



k-*0 



\ [tr[Cl] + (1 ± 1) tr[AC2]] =cl + ^^iT[AC\ (123) 



d(A; 



2^2 



tr[Z2] ± tr[Z 



21 



Ttr[Co^]tr[y^] 



1 tr[ACY + 2tr[Cg]tr[AC^] ltr[AC2]2 



21 



2 ti[Q?] ' 2 ti[n 

'^l^^l^'^l^^l±t«(-t.M.^l!^|. (124) 



Recall (see section 4.4) that the first of the physical conditions CI is equivalent 
to the statement that the 2x2 matrix Cq has two identical eigenvalues. But since 
Cq is symmetric this then implies Cq = Cq I, in which case the second condition is 
automatically satisfied. This also leads to the useful facts 

C/ab = =^ A = -2VPA0 PBO Uab; (125) 

2 UaaPao Ubb Pbo /ioc\ 

Co = = . (126) 

niA ruB 

Now that we have the fine tuning condition for the laser coupling we can compute the 
magnitude of the effective mass of the massive phonon and determine the values of the 
Lorentz violation coefficients. In particular wc shall start checking that this regime 
allows for a real positive effective mass as needed for a suitable analogue model of 
quantum gravity phenomenology. 
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6.1. Effective mass: 

Remember that the definition of mn reads 

m]i = hWjj/c^ (127) 
Using equation (125) and equation (126) we can rewrite Cq in the following form 

Co = [mBPAoUAA + ruAPBoUBB + UAB{pAomA + p Bom b)] / {Im Am b) ■ (128) 
Similarly equation (125) and equation (126) when inserted in equation (78) give 
2 4C/As(pAomB + pbo"^a)co 

^ ■ u^yj 

We can now estimate mn by simply inserting the above expressions in equation (127) 
so that 

^2 ^UAsipApmA + pBomB)mAmB , . 

[rriBPAoUAA + ttlaPboUbb + UAB{pAomA + PBQmB)] ' 
This formula is still a little clumsy but a great deal can be understood by doing the 
physically reasonable approximation tua ^ tub = m and Pa ^ Pb- In fact in this case 
one obtains 

2^2 8C/ab , . 

" " [f/.. + 2C/.. + C/BBl- 

This formula now shows clearly that, as long as the mixing term Uab is small compared 
to the "direct" scattering Uaa + Ubb, the mass of the heavy phonon will be "small" 
compared to the mass of the atoms. Though experimental reahzability of the system is 
not the primary focus of the ciurcnt article, we point out that there is no obstruction in 
principle to tuning a 2-BEC system into a regime where \Uab\ ^ \Uaa + Ubb\- For the 
purposes of this paper it is sufficient that a small effective phonon mass (small compared 
to the atomic masses which set the analogue quantum gravity scale) is obtainable for 
some arrangement of the microscopic parameters. We can now look separately at the 
coefficients of the quadratic and quartic Lorentz violations and then compare their 
relative strength in order to see if a situation like that envisaged by discussions of the 
naturalness problem is actually realized. 

6.2. Coefficient of the quadratic deviation: 

One can easily see from (123) that the ^2 coefficients for this case take the form 

zu2,i =0; (132) 
W2,ii = tr[AC'] = tr[Xi/'AX^/2] ^ ^^^^A] 
1 A f mAPAo + 'mBPBo\ 



2mAmB V y/PAoPBO J ' 

So if we insert the fine tuning condition for A, equation (125), we get 



(133) 



^2,11 UAB{rnAPAO + mBPBo) 
^2,/7 — — — 2 ■ \^^^) 

c^ mATnBC^ 
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^2,// = 774 . ^ , . ^ (135) 



Remarkably we can now cast this coefficient in a much more suggestive form by 
expressing the coupling Uab in terms of the mass of the massive quasi-particle mj^. 
In order to do this we start from equation (129) and note that it enables us to express 
Uab in (134) in terms of uj]j, thereby obtaining 

PAQiriA + pBomB u;jj 
4c^ PAoms + PBoruA mAms 
Now it is easy to see that 

PMmA + PBoms ^ 
PAoms + PBorriA 

and that this factor is identically unity if either niA = or pAo = Pbo- All together 
we are left with 

where 77 is a dimensionless coefficient of order one. 

The product in the denominator of the above expression can be interpreted as the 
geometric mean of the fundamental bosons masses uia and ms- These are mass scales 
associated with the microphysics of the condensate — in analogy with our experience 
with a 1-BEC system where the "quantum gravity scale" K in equation (5) is set by 
the mass of the BEC atoms. It is then natural to define an analogue of the scale of the 
breakdown of Lorentz invariance as M^s = -^mATnB- (Indeed this "analogue Lorentz 
breaking scale" will typically do double duty as an "analogue Planck mass".) 

Using this physical insight it should be clear that equation (137) effectively says 

*."-(S)^- (^^«' 

which, given that mj = 0, we are naturally lead to generalize to 

f f^i'xS^ / mass scale of quasiparticle \ ^ / or.\ 

'^^'^ = \ effective Planck scale j ' = ^'^^^ 

The above relation is exactly the sort of dimensionless ratio {p/My that has been 
very often conjectured in the literature on quantum gravity phenomenology in order to 
explain the strong observational constraints on Lorentz violations at the lowest orders. 
(See discussion in the introduction.) Does this now imply that this particular regime 
of our 2-BEC system will also show an analogue version of the naturalness problem? 
In order to answer this question we need to find the dimensionless coefficient for the 
quartic deviations, r/4, and check if it will or won't itself be suppressed by some power 
of the small ratio mu/Mes. 
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Let us now consider the coefficients of the quartic term presented in equation (124). For 
the various terms appearing in (124) we get 



tr[Z=^] = 2tr[M] = 
tr[AC2] _ i^[XK] -- 

tr[F2] = 2tr[XS-^] 



2 



m. 



m. 



(140) 



A rUAPAO + rUBPEO rr ^aPao + mspBO ... 

-77 , = 'Jab ; (141) 

Pao^ii.iUaa + Puo'iibUbb 



(142) 



2 PAomAPBOrriBUAAUBB 

where in the last expression we have used the fact that in the current scenario Uab — 0. 
Now by definition 



.(dA;2)2 

is the dimensionless coefficient in front of the k^. So 



fe=0 



?74 



h? 



± tr[Co^] 



tr[F 



21 ^ tr[AC2 



tr[fi2 



tr[Z2] ± tr[Z2 



2tr[Co2] 



± 



tr[F2] ^ tr[AC2] 



Whence 



M2 r 



+ 



tr[y^] 

2 

tr[AC^; 
tr[n2] 



+ 



2 ■ tr[l]2] 

tr[AC^] 
tr[Q2] 



Let us compute the two relevant terms separately: 



tr[Z^] 

tr[AC^] 
tr[02] 



4r2 



mAnii 



m 



A 



m 



rriAmB 



(143) 

(144) 
(145) 

(146) 
(147) 

(148) 



4M: 



PAprriAUAA + pBorriBUBB 
\_PaoPboUaaUbb [Uaa + Ubb j 

h? m\UAA + m^ilBB 



4M2g cl 



rnAruB {Uaa + Ubb^ 



(149) 



where we have used pxoUxx — n^'xcl for X — A, B as in equation (126). Note that the 
quantity in square brackets in the last line is dimensionless. So in the end: 



1 
4 



vriAmB 



m\UAA + m%UBB 
mATTiB (Uaa + Ubb 



m\UBB + m%UAA 
mArriB [Uaa + Ubb 



m\UAA + m^UsB 



mATUB 



(Uaa + Ubb^ 



; (150) 



(151) 
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Note: In the special case = we recover identical quartic deviations 774^/ = 
?74,// — 1/4, indicating in this special situation a "universal" deviation from Lorentz 
invariance. Indeed we also obtain 774^7 = 774^77 if we demand Uaa = Ubb, even without 
fixing niA = tub- 



6.4- Avoidance of the naturalness problem: 



We can now ask ourselves if there is, or is not, a naturalness problem present in our 
system. Are the dimensionless coefficients 774^7 /77 suppressed below their naive values by 
some small ratio involving M^s = ^JrrT^mB ? Or are these ratios unsupressed? Indeed 
at first sight it might seem that further supression is the case, since the square of the 
"effective Planck scale" seems to appear in the denominator of both the coefficients (150) 
and (151). However, the squares of the atomic masses also appear in the numerator, 
rendering both coefficients of order unity. 

It is perhaps easier to see this once the dependence of (150) and (151) on the 
effective coupling U is removed. We again use the substitution Uxx — "itt^xCq/ Pxo for 
X = A,B, so obtaining: 



rriApBO + TUB P AO 



m\pBo + ■m%pAQ 



(152) 



(153) 



jUAmB {itiaPbo + ruBpAo). 
From these expressions is clear that the 7/4 7 /77 coefficients are actually of order unity. 

That is, if our system is set up so that 777.77 ^ n^A/B — which we have seen in this 
scenario is equivalent to requiring Uab ^ Uaa/bb — no naturalness problem arises as 
for p > nijj Cq the higher-order, energy-dependent Lorentz-violating terms {n > 4) will 
indeed dominate over the quadratic Lorentz-violating term. 



7. Summary and Discussion 

Analogue models of gravity have a manifold role which can be summarized in three 
main points: (1) They can reproduce in a laboratory what we believe are the most 
important features of QFT in curved spacetimes; (2) they can give us the possibility 
of understanding the phenomenology of condensed matter systems via the body of 
knowledge developed in semiclassical gravity; finally (3) they can be used as test 
fields and inspiration for new ideas about the nature and consequences of an emergent 
spacetime [48, 49, 50, 51]. 

In this paper we have followed the last path by studying an analogue system 
which allows us to test the conjectures that lie at the base of most of the extant 
literature on quantum gravity phenomenology — by building an analogue spacetime 
exhibiting Planck-suppressed Lorentz violations. This analogue model, arising from 
a coupled 2-BEC system (previously studied in [39, 40]), reveals itself as an ideal 
system for reproducing the salient features of the most common ansatze for quantum 
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gravity phenomenology. Excitations in a coupled 2-BEC system result in the analogue 
kinematics for a massive and a massless scalar field. For a fine tuning in the 
hydrodynamic limit — sufficient to describe low-energy excitations — we recovered 
perfect Lorentz invariance. 

To describe highly energetic modes we modified the theory by including the 
quantum potential term. This is a quantum correction to the classical mean-field, 
which is energy-dependent and therefore is no longer negligible at high energy. We 
developed several mathematical tools to analyze the dispersion relations arising in this 
system. We Taylor-expanded the system around A; = 0, and calculated the coefficients 
for the quadratic and quartic order terms. We considered first the hydrodynamical 
approximation (phonons of long wavelength) as this case in a 1-BEC system leads 
unequivocally to a special relativistic kinematics. In the 2-BEC system we found that 
only for some specific fine-tuned values of the laser coupling A is a single background 
relativistic kinematics ( "mono-metricity" ) for the two phonon modes recovered. This is 
as expected as it is well known that for complex systems mono-metricity is not the only 
outcome for the propagation of linearized perturbations [49, 50, 51]. 

We then relaxed the approximation of long wavelengths for the 2-BEC system 
excitations. This allows to consider those Lorentz violations which are really due to the 
UV physics of the condensate, and which are indeed the source of the standard quartic 
Lorentz- violating term in the Bogolubov dispersion relation (5). By using an eikonal 
approximation we then studied the situation in which we again enforced the complete 
suppression of all the Lorcntz-violating terms which are not explicitly due to the UV 
physics, (of course this implies a fine tuning that is identical to that in the hydrodynamic 
approximation) . 

In the analogue spacetime so obtained we found that the issue of universality is 
fundamentally related to the complexity of the underlying microscopic system. As long 
as we keep the two atomic masses tua and ms distinct we generically have distinct 774 
coefficients (and the 772 coefficients are unequal even the limit tra — ruB)- However we 
can easily recover identical 774 coefficients, for instance, as soon as we impose identical 
microphysics for the two BEC systems we couple. Even more interestingly we saw that 
due to the presence of the interactions between the two BEC components, the quantum- 
potential-dependent Lorentz violations not only induce terms at order k"^, but also at 
order k'^ (in close analogy with what was predicted in [29] for a generic EFT with higher 
order Lorentz violations). Remarkably, we find that the 772 coefficients are exactly of 
the form envisaged (within the context of standard quantum gravity phenomenology) 
in order to be subdominant with respect to the higher-order ones. This implies that 
our 2-BEC analogue spacetime is an explicit example where the typical dispersion 
relations used in quantum gravity phenomenology studies are reproduced, and that 
the naturalness problem does not arise. 

Let us now comment briefiy about the nature of this result. First, it is important to 
stress that we did not merely perform a "tree level" calculation in the quasi-particle EFT. 
The dispersion relations we obtained were computed directly from the true physical 
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microscopic field theory describing the 2-BEC atomic system, and as such they already 
consistently take into account all the corrections due to higher loops in the quasi- 
particle EFT. In this sense the modified dispersion relations we found are those one 
would expect to observe if an actual experiment is set up, much as the form and 
the coefficients of the Bogoliubov dispersion relation are experimentally confirmed by 
present-day experiments [52]. Second, we stress that the avoidance of the naturalness 
problem is not related directly to the fact that we tuned our system so to reproduce 
special relativistic dispersion relations in the hydrodynamic limit. In fact our conditions 
for recovering SR at low energies do not a priori fix the the 772 coefficient, as its strength 
after the "fine tuning" could still be large (even of order one) if the typical mass scale 
of the massive phonon is not well below the atomic mass scale. Indeed the smallness of 
the r]2 coefficient is in this sense directly related to the mechanism providing a mass to 
one of the two phonons as we shall discuss at length below. 

The question we now want to address is why our model escapes the naive predictions 
of large Lorentz violations at low energies? There is a nice interpretation of this result 
in terms of "emergent symmetry" . 

We have seen that a non-zero A simultaneously produces a non-zero mass term for 
one of the phonons, and a corresponding non-zero LIV at order k'^ (single BEC systems 
have only LIV as described by the Bogoliubov dispersion relation). Let us now 
imagine driving A ^ but keeping the conditions CI and C2 valid at each stage (this 
requires Uab as well). In this case one gets an EFT theory which at low energies 
describes two non-interacting phonons propagating on a common background (in fact 
?72 — >^ and c/ = c// = Cq). This system posses an SO {2) symmetry corresponding to 
the invariance under rigid rotations of the doublet formed by the two phonon fields. 
Recovery of SR at low energies could then be seen as a by-product of imposing this 
symmetry on the system of two massless phonons. Hence non zero laser coupling 
A corresponds to a soft breaking of the S0(2) symmetry and of the corresponding 
Lorentz invariance at low energies. Such violation is then expected to be determined 
(as usual in EFT) by the ratio of the scale of the symmetry breaking m// and that of 
the scale originating the LIV in first place Mliv. We stress that the SO (2) symmetry 
is an emergent symmetry as it is not preserved beyond the hydrodynamic limit: the r/4 
coefficients are in general different if niA 7^ so 5*0(2) it is generically broken at high 
energies. However this is enough for the protection of the Zow;esi-order LIV operators. 

Hence this analogue model seems to be telling us that in solving the naturalness 
problem a possible mechanism could be that of an EFT which in the low-energy limit 
exhibits an accidental/emergent symmetry, instead of a fundamental one. Note that 
this is an interesting and original suggestion which goes in the opposite direction with 
respect to other attempts at solving the naturalness problem by looking at symmetries 
which are supposed to be exact in the high energy regime (see, e.g., works exploring 
the role of SUSY, such as [30], and the related problems with the low-energy symmetry 
breaking) . The lesson that can instead be drawn in this case is that emergent symmetries 
could be sufficient to minimize the amount of Lorentz violation in the lowest dimension 
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operators of the EFT. 

It is intriguing to think that an interpretation of SUSY as an accidental symmetry 
has indeed been considered in recent times [53], and that this is done at the cost of 
renouncing attempts to solve the hierarchy problem in the standard way. It might 
be that in this sense the smallness of the particle physics mass scales with respect 
to the Planck scale could be directly related to smallness of Lorentz violations in 
renormalizable operators of the low-energy effective field theory we live in. We hope to 
further investigate these issues in a future work. 
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Appendix A. Some matrix identities 

To simplify the flow of argument in the body of the paper, here we collect a few basic 
results on 2 X 2 matrices that are used in our analysis. 

Appendix A.l. Determinants 
Theorem: For any two 2x2 matrix A: 

det(A) = ^{tr[A]^-tr[A^]}. (A.l) 
This is best proved by simply noting 

det(A) = A1A2 = \ [(Ai + A2)' - (A^ + A^)] = 1 {tr[Af - tr[A^]] . (A.2) 

If we now define 2x2 "trace reversal" (in a manner reminiscent of standard GR) by 

A = A- tT[A] I; A = A; (A.3) 
then this looks even simpler 

det(A) = ~tr[A A] = det(i). (A.4) 

A simple implication is now: 

Theorem: For any two 2x2 matrices A and B: 

det(A + A 5) = det(yl) + A {tr[A]tr[5] - ix{A B\} + A^ det(5). (A.5) 
which we can also write as 

det(A + A S) = det(A) - A tr[>l B\ + A^ det(S). (A.6) 
Note that tr[A B\ = tr[A B\. 
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By iterating this theorem twice, we can easily see that: 
Theorem: For any three 2x2 matrices A, and C: 

det[A + A S + C] = det[A] - A ir{A B} + A^ [det[B] - tr{A C}] 

- \^ tT{B C} + X'^ det[C]. (A.7) 

Appendix A. 2. Hamilton-Cayley theorems 

Theorem: For any two 2x2 matrix A: 

iv[A] \-A A 



det[A] det[A]' 
Theorem: For any two 2x2 matrix A: 



(A.8) 



^tr[A] ±2Vdet^ j 
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